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G\ ; Abstract 

The coupling Qb^b-k is related to the form factor at zero momentum of the axial 
Ph! current between B*- and B-states. This form factor is evaluated on the lattice 

using static heavy quarks and light quark propagators determined by a stochastic 
inversion of the fermionic bilinear. The gB*Bn coupling is related to the coupling 
g between heavy mesons and low-momentum pions in the effective heavy meson 
chiral lagrangian. The coupling of the effective theory can therefore be computed 
by numerical simulations. We find the value g = 0.42 (4) (8). Besides its theo- 
retical interest, the phenomeno logical implications of such a determination are 
discussed. 
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1 Introduction 



The precise determination of the Cabibbo-Kobayashi-Maskawa (CKM) matrix will provide 
a stringent consistency test of the Standard Model (SM), together with new handles to 
understand CP violation and search for hints of new physics. As far as the heavy flavour 
sector is concerned, a large amount of experimental data is expected from I?-factories and 
CLEO in the near future. Non-perturbative QCD effects are the main source of systematic 
error in the extraction of fundamental parameters from experimental data: reliable results 
depend on some way to tame the large-distance dynamics. 

Lattice simulations provide a powerful tool to investigate non-perturbative dynamics 
from first principles, but systematic errors are introduced by the finite lattice spacing and 
the restricted range of quark masses which may be simulated. A different approach is based 
on exploiting exact or approximate symmetries of the theory to write effective lagrangians 
describing the large distance behaviour in terms of effective meson fields. The chiral sym- 
metry of strong interactions in the limit where m q — > constrains the terms appearing in 
the chiral lagrangian. The couplings in this lagrangian are phenomenological inputs, which 
need to be taken either from experimental data or from some other source. On the other 
hand, when the quark masses tend towards infinity, heavy quark effective theory (HQET) 
has proved to be a powerful tool to study heavy-flavour physics. However, HQET is less 
powerful when applied to heavy-to-light transitions, like B — > n, where the normalisation of 
matrix elements is not fixed by the symmetry. A combination of these two symmetries has 
been used in recent years to develop the heavy meson chiral lagrangian (HM%), describing 
the interactions of low-momentum pions with mesons containing a single heavy quark 
(for reviews see 0, [|). 

The definitions and notations used throughout this paper are conveniently introduced by 
a succint description of the building blocks of this effective theory. 

For the heavy degrees of freedom, heavy quark symmetry predicts that the wave function 
of a heavy meson is independent of the flavour and spin of the heavy quark, leading to a 
covariant representation of heavy mesonic states j3J. The angular momentum j and the 
parity P of the light degrees of freedom determine a degenerate doublet of states with spin- 
parity J p = (j ± l/2) p . The pseudoscalar and vector mesons (e.g. B and B*) correspond 
to the j = 1/2 case and are described by a 4 x 4 Dirac matrix H with two spinor indices, 
one for the heavy quark and one for the light degrees of freedom. In terms of the effective 
meson fields: 



where v is the velocity of the meson and B and B* are the annihilation operators for particles 
containing a b quark in the initial state. These meson fields are used in the effective lagrangian 
to describe the heavy mesons. The light mesons are treated as an octet of pseudo-Goldstone 
bosons according to the usual chiral lagrangian formalism. At low-momentum the strong 
interactions of B and B* mesons with light pseudoscalars are described by the couplings in 
the effective lagrangian; the lowest order interaction is given by [l], || : 





(1) 



= gTr(H a H b A b ;^ 7 5 ) 



(2) 



where 




(3) 
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with £ = exp(iM//). .M is a 3 x 3 matrix of 7r, 77 and fields 



(4) 



and the trace is over Dirac indices. At tree level / can be set equal to f n (the definition of 
the decay constant used here sets f n = 132 MeV). The roman indices a and b denote light 
quark flavour and repeated indices are summed over 1, 2, 3. The expansion of A in terms of 
pion fields begins with a linear term, 



A. = —A M + ■■■ ( 5 ) 



i 
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The coupling g in Eq. || can easily be related to the B*Bn coupling defined as [|| f|: 

(B°(p)n+(q)\B* + (p')) = -g B *BM 2 )q^{p') (2vr) 4 8(p' - p - q) (6) 

where rf is the polarisation vector of the B* and the physical states are relativistically 
normalised: 

(B(p)\B(p')) = 2p°(27r) 3 5 (3) (p-p') (7) 

The physical coupling gB*Bn is given by the value of the above form factor for an on-shell 
pion: 

9b*Btt = Jim QB*Bir{q 2 ) (8) 
At tree level in the heavy meson chiral lagrangian, the above matrix element is: 

(B (p)7r+(q)\B*+(p')) = g q^(p') (2vr) 4 5{p' - p - q) (9) 

which therefore yields: 

2rriB 

9b*b-k = —^9 (10) 

As a result, g and gB*B-w are considered as equivalent throughout this paper. The above 
relation can be extended to take into account higher-order terms in the HMx lagrangian J7L 
II- 

Starting from Eq. ^ and performing an LSZ reduction of the pion field, the coupling g is 
related to the form factor at zero momentum-transfer of the axial current between hadronic 
states. Such a relation is the analog, in the Bit system, of the Goldberger-Treiman relation, 
relating the nucleon electromagnetic form factor to the nucleon-nucleon-pion coupling. An 
important consequence of the Goldberger-Treiman relation, for our purposes, is that it allows 
a numerical evaluation of the gB*Bn coupling, as the form factors of the axial current can 
be evaluated by a lattice simulation. The details of the pion reduction are presented in 
Section 0. 

The interest of such a computation is two-fold. From a theoretical point-of-view, it is 
interesting per se to be able to fix, from lattice QCD, the coupling appearing in the heavy 
meson chiral lagrangian. On the other hand, it is important to stress that this determination 
also has phenomenological motivations. Assuming vector dominance in the B — > nlu decay, 
the coupling gg* s n fixes the normalisation of the form factors used to parametrise the matrix 
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element of the weak vector current, = wy^b, between hadronic states. Defining the form 
factors by: 

(n+(p>)\V»\B(p)) = h{q 2 ){p + p' •-^^-qr + 

= /+(<? 2 )(p+?/r+/-(<?v (ii) 

where q = p — p' is the transferred momentum, the contribution from the B* channel is easily 
evaluated: 

W,)= *WR (12) 

The normalisation of the form factor therefore depends on the B*B7r coupling and the decay 
constant of the vector meson, defined as: 

(0|V|B*(p)> =i^V(p) (13) 

JB* 

Heavy quark symmetry and chiral symmetry justify this pole form for fx when q 2 is close to 
9ma X = ( m B — m ir) 2 [0; B B B OH W ■ F° r <? 2 f ar fr° m ?m H ) the pole form may be taken as 



a phenomenological ansatz. However, we note that the functional dependence of the form 
factor in Eq. [12| cannot be simultaneously consistent with heavy quark symmetry at large 
q 2 , which demands that A^max) ~ m !/ 2 an d the light-cone sum rule scaling relation at 
q 2 = 0, which states fi{q 2 =0) ~ m B 3 ^ 2 fll2f . Nonetheless, by fitting lattice results, which 
are available in the high q 2 region where the pole form is justified, we can determine the 
parameters in Eq. 

It is interesting to remark that the interplay of the effective lagrangian approach and 
lattice simulations provides another determination of the form factors for the heavy-to-light 
B decays and therefore sheds further light on the theoretical determination of the non- 
perturbative effects mentioned at the beginning. This result can be compared with direct 
computations of the same quantities obtained by fitting lattice data [0, using unitarity 



bounds |14| and sum rules [g, [H| . 

In the work described here, the matrix element of the light quark axial current between 
the heavy mesons is computed in a quenched lattice simulation of QCD in the static heavy 
quark limit, using stochastic methods to compute the desired light quark propagators, as 
described in Section |3|. The discussion of systematic errors is an important issue in any 
lattice calculation and plays a crucial part in estimating the error on the final result. In this 
respect, it is important to stress here that we are presenting an exploratory study. Our main 
concern is therefore to test the possibility of extracting the coupling defined above, rather 
than presenting its best lattice determination. Such a task would require a more extensive 
simulation and is left for further studies. 

Renormalisation constants are needed in order to connect lattice results with continuum 
physical observables. Those relevant for the action and the quantities considered in this 
paper are summarised in Sect. ^. 

The best estimate we obtain for g is 

<7 = 0.42(4)(8) (14) 
The phenomenological implications of this result are discussed in Sec. [5[ 



3 



2 Pion reduction 

An LSZ reduction of the pion in the definition of Qb* b-k > Eq. || yields: 

(B°(p)n + (q)\B* + (p + q)) = i{ml - g 2 ) / e^(B(pMx)\B*(p + q)) (15) 
Using the PCAC relation: 

<x) = 4f^(i) (16) 

mi Jtt 

where is, as usual, the QCD axial current, Eq. [15| becomes: 

(B\p)ir + (q)\B* + (p + q)) = q^™tL^ f e ** (B(p)\A,(x)\B*(p + q)) (17) 

Jtt Jx 

The matrix element of the axial current is parametrised in terms of three form factors: 



(B\p)\A li mB* + {p + q)) = V,F l (q 2 ) + ( V -q)(2p + q) ll F 2 (q 2 ) 

+ (v-q)q^ 3 ( q 2 ) (18) 

yielding for the B*Bn coupling: 

9b*bM 2 ) = -T 7 ^- \ F ^ 2 ) + {<• ~ m %) F ^ 2 ) + ^fa 2 )] ( 19 ) 

fix m jr ' J 

In the static limit in which our simulation is performed, the B and B* mesons are degenerate, 
so that the form factor F 2 can be discarded. 

Analytical continuation of Eq. [19| towards the soft-pion limit (q 2 —>■ 0), leads to: 

9b*bM = - T Fi{0) (20) 

Jtt 

It is commonly assumed, when deriving the Goldberger-Treiman relation, that Qb^Btt is a 
smooth function of q 2 , and, therefore, that the physical coupling can be approximated by: 

9B*Btt = gB*BTT( m l) ~ 9B*Btt(0) (21) 

The above equation explicitly shows that, in the soft-pion limit, the B*Bir coupling is related 
to the form factor of the axial current between B and B* states. If one were working in the 
chiral limit from the very beginning, the same Goldberger-Treiman relation, Eq. would 
be obtained from the conservation of the axial current. 

The relation between gs*BTT and g mentioned earlier can be rederived by comparing the 
matrix element of the Noether current associated to chiral symmetry both in HMx and QCD. 
In the chiral limit, the Noether currents associated with chiral symmetry, in QCD and in 
HMx, are respectively 

3c&* = ?W*f (22) 

and 

Jhm x „ = fnd,M ab - 2g (B; a ^ B b + B a *B; b ) + ■ ■ ■ (23) 

where the ellipsis denotes terms with more than one pion or terms containing both heavy- 
mesons and pions. 
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Figure 1 Tree-level diagrams needed to compute the matrix element of the axial current in HM^. 
The dot (•) represents the g- vertex in the HM^ lagrangian; the squares are insertions of the axial 
current in Eq. [23|. 



The form factors of the axial current in QCD are related to the coupling of the heavy 
meson chiral lagrangian by matching the two theories at tree level. The diagrams needed 
at tree level to evaluate the matrix element in Eq. |T^ for the HMx current are depicted in 
Fig. |I[ A straightforward computation leads to: 

F 1 (q 2 ) = -2gm B 

2g 



W) 



Q 



-m B 



which, together with the Golberger-Treiman relation, reproduces Eq. [10]. One may also work 
away from the chiral and heavy quark limits and include corrections for finite mass pions 
and heavy quarks in this result. 

The determination of g is therefore reduced to the computation of the matrix element of 
the light-light axial vector current between hadronic states. Such an evaluation can be per- 
formed using three-point correlation functions on the lattice. The details of the calculation 
are reported in the next section. 



3 Lattice results 



3.1 Stochastic propagators 

The numerical analysis is carried out on 20 quenched gauge configurations, generated on a 
12 3 x 24 lattice at — 5.7, corresponding to a -1 = 1.10 GeV. The heavy quark propagators 
are evaluated in the static limit |17|]. Stochastic propagators fjl| are used to invert the 
fermionic matrix for the light quarks. They can be used in place of light quark propagators 
calculated with the usual deterministic algorithm. The stochastic inversion is based on the 
relation: 



i.i 



Ml 



1 

~Z 



V<f) (M ifc fc )*& exp (-#(M+M) y 0. 



(24) 



where, in our case, M is the tadpole improved SW fermionic operator and the indices i, j, k 
represent simultaneously the space-time coordinates, the spinor and colour indices. Two 
values of k are considered, K\ = 0.13843 and k 2 = 0.14077, with c sw = 1.57. The lighter 
value k>2 corresponds to a bare mass of the light quark around the strange mass. The chiral 
limit corresponds to k c = 0.14351 | 20| . For every gauge configuration, an ensemble of 24 
independent fields 0j is generated with gaussian probability: 



P[0] = lexp (-0*(MtM) 



(25) 
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All light propagators are computed as averages over the pseudo-fermionic samples: 

Sij = \ or (26) 
{ ((750*),- (M075),) 

where the two expressions are related by 5 = 756^75. Moreover, the maximal variance 
reduction method is applied in order to minimise the statistical noise [jnj. Maximal variance 
reduction involves dividing the lattice into two boxes (0 < t < T/2 and T/2 < t < T) and 
solving the equation of motion numerically within each box, keeping the spinor field (p on 
the boundary fixed. According to the maximal reduction method, the fields which enter 
the correlation functions must be either the original fields <fi or solutions of the equation of 
motion in disconnected regions. The stochastic propagator is therefore defined from each 
point in one box to every point in the other box or on the boundary. For this reason, when 
computing a three-point correlation function 

< 1 , a;) O(t , y) «/ f (^2 , ^) 1 0> (27) 

one operator — O in the present work — is forced to be on the boundary (t = or T/2) and 
the other two must be in different boxes, while the spatial coordinates are not constrained. 
If j is a point of the boundary, not all the terms in (M(j))j lie on the boundary because the 
operator M involves first neighbours in all directions. Hence, whenever a propagator Sij is 
needed with one of the points on the boundary, we use whichever of the two expressions in 
Eq. |2B| has the spinor M<fr computed away from the boundary. 

In smearing the hadronic interpolating operators, spatial fuzzed links are used. Following 
the prescription in |TB|, to which the interested reader should refer for details, the fuzzed 
links are defined iteratively as: 

U new = V (fU old + Ubend,^ (28) 

where V is a projector over SU (3) and £/b en d,i are the staples attached to the link in the 
spatial directions. We take / = 2.5 and use two iterations with fuzzed links of length one. 



The smeared fermionic fields are defined following ET] . 



3.2 Lattice computation 

In order to extract g, the three-point correlation function C3 and the two-point correlation 
function C2 for local (L) and fuzzed (F) sources are computed. The FF three-point function 
is defined as: 

CiU^MM = ^jd'y (0|Jf (y.-tO^Cx + y.oy* t(y,t 2 )|0> (29) 

where and J B are fuzzed operators with the quantum numbers corresponding respec- 
tively to the B* and B states. Analogous definitions hold for the FL and LL cases. For 
time separations that are large enough to isolate the lowest energy states, the three-point 
function is related to the axial current matrix element: 

cii^-MM) - (our \k) { -^^ (b\j b me- MB{t ^ 

= ^^( ) gM / e -^(wi 2 ) (30) 

2/71 £J 



6 



where r is the polarisation label of the vector particle, and the sum over polarisations is 
omitted. In the static limit considered in this paper, the slope of the exponential time-decay 
is not the physical mass of the mesons; it can be interpreted as a binding energy. Moreover, 
the two-point functions for the vector and pseudoscalar particles are degenerate, leading to 
the same binding energies and Z factors for both the B and B*. In order to avoid confusion, 
the physical mass is denoted by and the binding energy by Mb- Z f is the overlap of the 
interpolating operator with the physical state, defined from the two-point functions: 

CZ F (t) = ^J v d 3 y(0\J B (y,0)J B \y,t)\0) 

-> (Z F ) 2 e- MBt (31) 

Integrating the three-point function over x gives the matrix element for zero momentum 
transfer. In this limit, the latter can be expressed in terms of the form factor Fi(q 2 ) in 
Eq. pi. 



The sum over polarisations in Eq. 30 yields: 



/• 



d^xC^M^^iZy^-^F^-"^' (32) 

The last equation shows that the three-point functions with /i ^ v and \i — v = should 
vanish. Therefore, only the correlators with \i — v— 1, 2, 3 are henceforth considered. 

Moreover, taking rotational symmetry into account, C F F (x; ti, £2) is expected to be a 
function of the distance r only, up to cut-off effects. The three-point functions measured on 
the lattice are averaged over equivalent x positions^. The desired matrix element is obtained 
from the ratio: 

E (r t) * (Z^V* ^ ( B *\Mr)\B) . , 

E,{r,t) - [Z ) ^ rC F F{t)c ? F{t) - ^ % (33) 

where the time-dependence cancels when the three-point function is divided by the product 
of two-point functions. 

The coefficients Z F and Z L are extracted from the exponential fit of the two point 
correlation functions C2. The data are reported in Fig. 0. As one can see from the plots, 
the two-point functions already exhibit a single-exponential behaviour at moderate time 
separations. The main sources of error in this computation are expected to stem from the 
determination of the three-point function and from the value of the light quark masses, 



1 The symbol 

E, 



indicates the average on all the spatial positions on the lattice compatible with the constraint \x\ = r. On a 
finite lattice V — L 3 , only some distances are allowed 



because each xi must be an integer between and L/2. For each allowed distance r, a given number of terms 
N r appears in the above sum, yielding a relative error on each point SE^r, t) ~ N r , e.g. 

N = 1; N t = 6; = 12; = 6; . . . ; = 1 
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Figure 2 Logarithmic plots of \t) and C2 L (t) for both values of /c. The quoted values refer 
to the reduced \ 2 ■ 
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M B a 


0.930(4) 


0.898(5) 


0.862(7) 


(Z F f 


4.75(10) 


4.57(10) 


4.37(15) 


(Z L f 


0.59(1) 


0.58(2) 


0.57(3) 


/| tatic (GeV) 


0.43(1) 


0.43(1) 


0.42(2) 



Table 1 Values for Z F , Z L , Mb and J^ atlc obtained from fitting the two-point functions. 

which are far from the chiral limit. Thus, a single exponential fit of the two-point functions 
turns out to be accurate enough for the scope of this study. The value of Z F is extracted 
from a direct fit of C FF , while Z L is obtained from C FF and C FL . The results of the fit 
are shown directly on the plot. It is worth remarking that, in all the channels considered, 
single-exponential fits yield reasonable values for the reduced \ 2 . 

The B meson decay constant, in the static approximation, can be extracted from: 

^static = ^static / J_ z L fl -3/2 ( 34 ) 
V m B 

where Z^ atlc is the renormalisation constant for the axial current in the static theory, which 



is discussed in the following section and defined in Eq. The aim here is not a precise 
determination of the pseudoscalar decay constant, f^ atlc . Rather, the result is presented to 
allow an estimate of the systematic errors in our computation of the HMx coupling, g. 

The results of the fits, together with the values for Jj* atlc are summarised in Tab. p]. The 
B meson binding energies obtained from the fits of C FF and C FL are consistent with each 
other. However, our determination appears to be slightly different from the one published 
in [IE . The discrepancy can be explained as a contribution from excited states, which is 
subtracted in |18j where a multi-exponential fit is performed. The 1-state fit yields a value 
of Z L approximately 20% higher than the one obtained from the 3-state fit. Hence the 
value obtained for the static B decay constant lies above other lattice calculations of this 



quantity jl6fl . It is striking that the actual number does not depend on the value of the 
hopping parameter k. However the variation, as the bare mass of the quark goes from m s 
towards the chiral limit, is also expected to be about 20% and could be obscured by the 
contamination from excited states. 

This is a first, exploratory, direct lattice determination of the coupling g, so the discrep- 
ancies noted above are perhaps expected and could easily arise from various lattice artefacts. 
Our value of (3 is far from the continuum limit, while our action and operators are not fully 
0(a) improved. We have not tried to optimise the smearing or fitting procedures. Our 
ensemble of gauge configurations and the collection of spinor configurations on each gauge 
sample are quite small. The calculation is also performed in the quenched approximation. 
All of these issues could be addressed in a further simulation, but here we will keep in mind 
that the uncertainties will propagate as systematic errors to our final result for the B*Bir 
coupling. 

The generation of stochastic propagators for the light quark and the static approximation 
for the heavy quark have proved to be useful tools in this investigation. However, neither is 
strictly necessary to calculate the axial current matrix element of interest. One could com- 
bine static heavy quarks with light quark propagators determined by standard deterministic 
methods. A full 0(a)-improved simulation with heavy quark masses around the charm mass 
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would allow one to go beyond the static approximation and study the dependence of g on 
the heavy mass. The freedom allowed in lattice calculations to tune quark masses would also 
allow the light quark mass dependence, noted above as strikingly absent, to be investigated 
in more detail. 

Returning to the results of the present simulation, the quantity E (x, t), which is expected 
to vanish, is measured as a further consistency check. The data reported in Figs. |3| and |] 
show a much smaller signal than the one obtained for E i {x 1 1). As t is increased, the fictitious 
peak at zero distance decreases, while the noise increases. 

Using rotational invariance, the average of the three spatial components of Ei(r,t) 

E(r, t) = ^(r, t) + E 2 (r, t) + E 3 (r, t)) (35) 

is measured. The results are reported in Figs. |] and |], for the two values of k used in this 
simulation. At fixed r, E(r,t) is expected to be independent of t. As this behaviour is 
confirmed by the data, the signal can be improved by averaging the values at different times, 
each weighted with its error, yielding a function of the spatial distance E(r), which needs to 
be integrated over the three-dimensional volume. The time-slices considered in the average 
are t = 4, 5, 6. Logarithmic plots of E{r) is displayed in Figs. || and || for both values of k, 
suggesting that the data are consistent with an exponential decay. 

To evaluate the volume integral, E(r) is fitted with a two-parameter exponential for each 
value of k, 

f(r) = Se~ r/r ° (36) 

The results of the fit and the values of the reduced x 2 are recorded on the figures. The 
coupling g is finally obtained by integrating analytically the fitted function: 

g L = - A7rr 2 E(r)dr = - / 47rr 2 f(r)dr (37) 
Jo Jo 

The superscript L indicates that these numbers are defined on the lattice using operators 
renormalised at the lattice energy scale a -1 = 1.10 GeV. As for the two-point functions, 
the value of g h does not appear to depend on the mass of quarks. It is not clear from this 
simulation, whether this is a genuine physical feature or, as explained earlier, an artefact of 
the lattice used here. Further studies should aim to clarify this point. 

The best estimate for g L at k = k c is thus obtained from a weighted average of the results 
at the two kappas used in our simulations: 

g L = 0.52(5)(10) (38) 

The first error is statistical; the second is a conservative estimate of the systematic error of 
20%, based on the systematic error observed in the estimate of 



4 Renormalisation constants 

Quantities evaluated on the lattice are connected to their continuum counterparts by a finite 
renormalisation. In order to extract physical information, the matrix element of a bilinear 
quark operator defined on the lattice, C L (a), has to be multiplied by the corresponding 
renormalisation constant 

0(fi) = Z (afi,g)O h (a) (39) 
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Figure 3 Plots of E°(r,t) and E(r,t) as functions of r for different values of t (= 3,4,5,6) for 
k = 0.13843. 
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Figure 4 Plots of E°(r,t) and E(r,t) as functions of r for different values of t (= 3,4,5,6) for 
k = 0.14077. 
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Figure 5 Plot of E(r) as a function of r for k = 0.13843. 
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Figure 6 Plot of E(r) as a function of r for k = 0.14077. 
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which in general depends on the renormalisation scale fi and the details of the lattice dis- 
cretization (a single continuum operator may also match onto a set of lattice operators). 
For partially conserved currents, the \i dependence disappears in the above definition, the 
associated anomalous dimension being equal to zero. This is the case for the QCD light-light 
axial current considered here. It is also the case for the heavy-light current in full QCD, 
but not for the static-light current onto which it matches in the effective theory: hence the 
renormalisation constant Z s J^ tlc in Eq. [JU] below, which converts the lattice matrix element 
to the physical /| tatlc , includes the effect of running between different scales in the effective 
theory. 

In our simulation, we use the standard gluon action and the tadpole-improved p2| 



Sheikholeslami-Wohlert [23] fermion action for the light quarks: 



a 4 



S = ^-—^^ v {x)5 xy (40) 



auQ 



(1 - -i^)U^(x)5 x+ ^ y + (1 + lp)UtWx-[i, v w{y) 



1 



where as usual F^ v (x) is a lattice definition of the gluon field strength tensor, k = kuq and 
we define uq from the average plaquette in infinite volume, Uq = (| Tr U a ). The redefinition 
of the parameters in the action is done in order to absorb the large renormalisation coming 
from lattice tadpole graphs [P^| . The standard definition of c sw , the coefficient of the u^ v F t 
term, in this case gives c sw = 1/uq. For our lattice, u = 0.86081 from the average plaquette 
so that c sw = 1.57. For the b quarks we use the standard static quark action. Here we 
summarise the required renormalisation constants for the light-light and static-light axial 
currents with the action and parameters defined above. 

The matrix element of the light-light axial current between an initial state % and a final 
state / in the continuum can be written as: 

(f\A,\i) = Z A (f\Al\i) = Z^^^iflA^i) (41) 

Uq 

The factor it , a measure of the average link variable, can be interpreted as a non-perturbative 
rescaling of the quark fields in the tadpole improvement prescription. The 1-loop perturba- 
tive renormalisation constant must then be redefined by dividing out the 1-loop expression, 
w o~ lo ° P ) corresponding to the chosen definition of uq, 

r^tadpolc ^1— loop / 1— loop //1o\ 

L A — L A l u l 4Z J 

The overall renormalisation constant is given by the whole factor 

Z A = u ^ adpole (43) 

The perturbative part reads: 

gr^ = i+f cx „ 1+ fg£ (Cx _ A) (44) 



where A = —tt 2 for our plaquette definition of u and (a = —13.8 |25|, |26], ^7 . 

We note here that the removal of tree- level 0(a) discretisation errors is achieved by 
combining the Sheikholeslami-Wohlert action for c sw = 1 with improved operators, found 
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by redefining or "rotating" the quark fields [^8]. In our simulation, the rotation has not 
been applied to the quarks, so the perturbative coefficient ( above is calculated using the 
improved action only. Moreover, in perturbation theory, c sw = 1 + 0(a), so we quote ( 
above calculated for c sw = 1, without rotated light quark fields. In [^J the light fermion 
bilinear operator renormalisations are given as functions of c sw and an arbitrary amount of 
field rotation: using our actual value of c sw = 1.57, with no field rotation, in those results 
would increase Z A by 5%. 

The mean-field improved perturbative expansion is performed in terms of a boosted 
coupling constant, a, which in our case is chosen as a = a /uQ, where a is the bare lattice 
coupling constant. We use a for a in Eq. [44|. For our lattice, with (3 = 5.7 and u = 0.86081 
from the average plaquette, this leads to: 

Z A = u Z A adpole = 0.806. (45) 

For the heavy- light axial current computed in the static limit, the renormalisation con- 
stant Z A tatlc is found by a three step procedure combining the matching between full QCD 
and the continuum static theory at a scale of order nib [0' the continuum static theory 



running from nib to a lattice scale q* [31, 32. 33 1 and finally the matching between the con 



tinuum static and the lattice static theories |34|, 35| . We extract from the calculation of 



Borrelli and Pittori [35[ the appropriate matching factor corresponding to our case of an 
improved action, but without field rotations on the light quarks. In the notation of [|35| the 
number we need is 

C F {5/4 - A 1Q15 - d 1 - f/2 + (e - e rcd )/2) = -19.36 (46) 

This appears in Eq. [19| below and contrasts with the values —27.2 for the Wilson action and 
—20.2 for the improved action with c sw = 1 and field rotations included. 



In the matching to the lattice we will adopt a Lepage-Mackenzie choice |22| of the scale 
q* at which to perform the matching. This scale is determined from the expectation value 
of ln(ga) 2 in the one-loop lattice perturbation theory integrals for the corrections to the 
renormalisation constant, including the perturbative tadpole improvement corrections for 
the chosen definition of u . For the improved, c sw = 1, action and plaquette definition of u , 
Gimenez and Reyes quote q*a = 2.29. We will use this value. 



We will also adopt a plaquette definition of the lattice coupling constant according to 
(for zero flavours): 

4:71 

-ln((Tr[/ D )/3) = — a y (3.41/a)(l - 1.19a v ) (47) 
3 

Once ay(3.41/a) is determined, we can use the equation for the running of ay, 

4tt 

^ = 2boHq/A v ) + b 1 In (2 ln(q/A v ))/b (48) 

to determine av (<?*). In the quenched theory, b = 11 and b\ = 102. We find aAy = 0.294 
and a v (q*) = 0.216. 

Since q* = 2.52 GeV is between the charm and b quark thresholds we perform the contin- 
uum running with four active flavours. From the Particle Data Group (PDG) ||37| , we take 
— 237^24 MeV using two-loop running, and find using continuity of the strong cou- 
pling at the b quark threshold. This threshold value is given by nib satisfying ^^(m;,) = nib. 
We take m b = 4.25 GeV, using the average of the range, 4.1-4.4 GeV, quoted by the PDG 
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The overall renormalisation constant is thus given by: 



^static / a(m b ) \ d r c a{m b ) a{g*)-a{m b ) J 
A \ a(q*) J I 4vr 471^ 

x ^(l + ^M[41n(g*a) - 19.36 - C F X/2]) (49) 

where c = 8/3, = 0.91 and d = —6/25. Using the inputs given above, we find 

^static = 7g ^ 

Independent variations of A^ 5 ^ and nib to the ends of the ranges quoted above and ±10% 
variation in a -1 change this value by 1.3% or less. Changing q* to 1/a or n/a reduces Z s J^ tlc 
by 13% or raises it by 1.3% respectively. 

We close this section by noting that both and Z^ atlc are evaluated in the chiral limit 
and so do not depend on the light quark mass used in the simulation. 



5 Phenomenological implications 

We now combine the lattice matrix element with its renormalisation factor to find the con- 
tinuum value for g. Since we observe no light quark mass dependence in g L and evaluated 
the renormalisation constant in the chiral limit, we multiply the lattice matrix element in 
Eq. |38| by the renormalisation constant in Eq. H5L to obtain the value: 



g = Z A g L = 0A2(4)(8) (51) 

for the coupling of the heavy mesons with the Goldstone bosons. 

The direct decay B* — > Bit is kinematically forbidden. However, the corresponding 
reaction occurs in the D system, where the coupling go* d-k is also related to g by an expression 
analogous to Eq. [10], although the 1 / uiq corrections are expected to be larger in the charm 
case. A recent analysis fl38 |, incorporating chiral symmetry breaking corrections plus 1/m. 



c 

corrections in the HMx lagrangian and fitting to the branching ratios for D*° —>■ D°n°, 
D* + — > D + n° and D* — > D s n° (together with radiative decay rates for the same D* mesons), 
gives 

o-|°- 27 ^® (52) 
9 ~ \0.76(3)(?) 

The two fold ambiguity can be resolved by imposing the experimental limit T(D* + ) < 
0.13 MeV [|%9|1 , which gives g < 0.52 to a good approximation |38| . 



Other estimates of g are derived using constituent quark models and sum rules. Starting 
from the non-relativistic result g = 1, quark models can be improved using more sophisticated 
assumptions to describe the quark dynamics inside the meson. Independent estimates are 
obtained from computing QCD correlation functions in the framework of sum rules. To give 
an idea of the range spanned by different determinations, some recent results are listed in 
Tab. g 

The best estimate from a global analysis of available results, quoted in the review @ , is 

g ~ 0.38 (53) 
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Reference g Reference g 

B 0.32(2) H 07T 



g 0.39 § 0.75-1.0 



0.44 (16) H, 0.4-0.7 

0.38 



a Combined sum rule + lattice results 

Table 2 Recent determinations of the coupling constant g 

with a total uncertainty of 20%; the latter being comparable with the estimated systematic 
error from our lattice computation. Not only is the result presented in this paper compatible 
with the previous estimates; the systematic error is also competitive when compared to the 
above results. 

The agreement with other previous estimates is pleasing, but the value of the coupling 
g can also be used to check the consistency of the heavy quark symmetry predictions in the 
soft pion limit for the lattice results. 

The form factor fx of the semileptonic decay B — > tcIu is predicted to have the expression 
given earlier in Eq. [12]. This behaviour for fi{q 2 ) is expected to be valid near zero recoil 
(q 2 — > <?ma X )' where the closest vector resonance dominates, even beyond the leading ap- 
proximation in 1 /ni}, m HQET JET]. One also obtains this behaviour from the heavy meson 



chiral lagrangian [1], [7], |8|, fj, valid for a low momentum pion. It is perhaps surprising that 
such a behaviour is found by sum rules to hold reasonably, at least for the D meson, also 
at q 2 -> f§. In this case the coupling g would fix the normalisation of the form factor 
fx . The extension of the validity of the vector pole dominance for general values of q 2 has 
no simple theoretical justification. It can be argued that the contributions from different 
resonances can lead to a single pole behaviour; however, in this case, the relation between g 
and /i(0) would be spoilt. 

Lattice data for the semi-leptonic B decay form factors have been fitted assuming a pole 
behaviour for f\ [O, yielding /i(0) = 0.44(3), to be compared with 

fi(0)\vun = ^f~9 (54) 

JttJB* 

from Eq. [12]. Our determination of g gives /i(0) = 0.50(5) (10), in good agreement with the 
direct fit. Such an agreement should not come as a surprise: the lattice data, after chiral 
extrapolation, turn out to be close to the end-point of the phase space kinematically allowed 
for B —>■ it decays (q 2 ~ qm ax )- The lattice normalisation of the form factor comes therefore 
from a fit in a region where vector dominance does hold. Hence, the above agreement should 
be seen as a consistency check of the two lattice computations. It is nevertheless important 
to see that the two results are not contradictory. 

Assuming a pole form for fx, an independent bound on the value of the residue is obtained 



by enforcing the theory to satisfy unitarity |45[| . The bound quoted in jE| is 

A(0)m|, < lOGeV 2 (55) 

which translates into 

/i(0)<0.4 (56) 

This, again, agrees reasonably with the determination obtained from VMD, using our value 
for g, and the one coming from the direct fit of the lattice form factors. 
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As the lattice determinations will improve in the future, the comparison of the three 
results summarised here could become an effective way to constrain the residue at the B* 
pole. 



6 Conclusions 

We have shown that the coupling gB*Bn, or equivalently the coupling g in the HM% la- 
grangian, may be evaluated on the lattice from the matrix element of the light quark axial 
current between B and B* states. The value of g enters many phenomenological quantities of 
interest calculated in heavy meson chiral perturbation theory, including the form factors for 
semileptonic B —>■ it decays mentioned here, as well as B* — ► 7r decays and other quantities 
such as ratios of heavy meson nonleptonic decay constants, heavy meson mass splittings and 
radiative decays. 

Even the relatively crude estimate obtained here should be interesting for heavy meson 
phenomenology, but we believe future more precise lattice computations would be valuable. 
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